N85-25890 


DIRECT TRANSIENT ANALYSIS OF A FUZE 
ASSEMBLY BY AXISYMMETRIC SOLID ELEMENTS 

Chung C. Dai and Jackson C. S. Yang* 

Advanced Technology & Research, Inc. 

Burtonsville. 55D 20866 

John Titus 

Harry Diamond Laboratory 
Adelphi. MD 20783 

ABSTRACT 

A fuze assembly. which consists of three major parts, i.e., nose, collar 
and sleeve, was designed to survive severe transverse impact giving a maximum 
base acceleration of 20.000 G. Experiments showed that hoop failure occurred 
in the collar after the impact. They also showed that by bonding the collar to 
the nose, the collar was able to survive the same impact. To fjrnd out the 
effectiveness of the bonding quantitatively, axisymmetric solid elements TRAP AX 
and TRIAAX were used in modelling the fuze and direct transient analysis was 
performed. The dynamic stresses in selected elements on the bonded and 
unbonded collars were compared. The peak hoop stresses in the unbonded collar 
were found to be up to three times higher than those in the bonded collar. 
NASTRAN results successfully explained the observed hoop failure in the 
unbonded collar* In addition, static and eigen value runs were performed as 
checks on the models prior to the transient runs. Their results are compared 
with experiments. Comments on the use of the MPCAX cards, the existence and 
contributors of the calculated first several nearly identical natural 
frequencies, are also addressed. 


INTRODUCTION 

The M732 E2 Artillery Proximity Fuze comprises an aluminum fuze body which 
houses a safety and arming module, detonator bulkhead, and power supply aft of 
a center bulkead. Forward of the bulkhead is the fuze processor electronics 
and loop oscillator. The oscillator is housed in a 30% glass reinforced 
modified Polyphenylene Oxide (PPO) nose cone which is secured to the fuze body 
with a 30Z glass reinforced PPO threaded retaining collar. See Figures 1 and 
2 . 


When a projectile is fired from a gun, the fuze is subjected to side loads 
approaching 20,000 G's. In the laboratory, a machine was developed to simulate 
this loading. The fuze is secured to a rail guided steel fixture wlich is 
struck from the side by a steel weight. The impact is tailored to give a 
20.000 G peak lateral acceleration with a 30 fps velocity change. 

While fuzes at ambient temperature experienced no r a u;ee vben subjected 
to the test, fuzes which were temperature conditioned at -40^, failed without 
exception. The major failure mode, named hoop failure, appeared to be hoop 
stress on the retaining collar causing a longitudinal crack from end to erd. 

See Figure 3. 
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It was discovered that by solvent bonding the retaining collar to the nose 
cone during assembly, the lov temperature failures could be avoided. Since the 
mechanism by which the solvent bond alleviated the problem was not well 
understood nor was the safety factor of the bond known, it was feared that 
future modifications to the fuse could not be evaluated with respect to the 
impact on the fuze joint integrity. Therefore. NASTRAN was used to model the 
forward section of the fuze in an attempt to better understand the behavior of 
the current design, ana better predict the behavior of future designs. 

NASTRAN MODEL 

The given fuze forward section consists of three major parts* i.e. f nose, 
collar and sleeve. As depicted in Figures 4a. 4b, and 4c, 105 TRAPAX and 9 
TRIAAX elements were used to model this fuze and its interior electronic parts. 
These elements are inter-connected at 180 different rings or circles. Static 
and eigen value analyses were performed prior to the final transient runs, 
numerical results from the static and eigen value runs were carefully examined 
and compared with the available experiment data. Constraint conditions and 
boundary conditons had been modified accordingly. We will discuss in more 
detail later. 


TEST MODEL I 

A cantilevered, hollow slender cylinder, as depicted in Figure 5a. was 
modelled by TRAPAX elements as a preliminary study. Stress analysis of this 
cylinder was performed. Axial pressure load or concentrated transverse tip 
loading was separately applied. Numerical results from these test cases 
provided basic understanding on TRAPAX's behavior and characteristics. Axial 
pressure loading condition was first considered, see Figure 5b« As expec:ed, 
NASTRAN results showed that the cylinder's cross sections were displaced 
uniformly and axial stress was developed uniformly across each cross section. 

Figure 5c shows the transverse tip loading condition. The cylinder was 
modelled by regular 8x2, 16x2. 32x1. 32x2. 64x1 meshes. Both the tip 
deflection and the bending stress at point A were checked for convergence with 
respect to mesh refinement* As noted :n Table 1. with the same number of 
elements. 64 x 1 mesh gave better resuls than 32 x 2 mesh. B?sed on this 
observation, we used only one layer across the thickness of the real model. 

TEST MODEL II 

A cantilevered, hollow slender cylinder of different size (Figure 6) was 
used in the remaining preliminary studies. It was modelled also by TRAPAX 
elements but with 4 x 1 or 5 x 1 mesh. This mode 1 served for testing NASTRAN 
rigid formats 3. 9, and 12. simulating the presecribed base acceleration, and 
also for testing NASTRAN checkpoint and restart features. 

Free vibration analysis was performed on the cantilevered hollow cylinder 
in Figure 6a. Different frequency ranges were specified requesting ths 
extraction of natural frequencies and modes. 
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Given V * 30 ft /sec after the impact and the aaximia acceleration of 
20,000g, the duration of the impulse depends on the assumption of its shape. 

For a sinwt shaped impulse, the duration is about 70 z 10 sec. For a sin*wt 
shaped impulse, the duration is about 90 x 10* sec. Impulses of both shapes 
mere tested on Model II. However, only the sin* cot shaped impulse was employed 
in the analysis of the real model because it simualtes the real impulse more 
closely. 


Several schemes to input the prescribed base acceleration a(t) were also 
carefully examined. First, a nearly rigid layer with huge mass M was attached 
to the bottom of the original cylinder. Dynamic force F(t) * Ma(t) was then 
applied to this layer in several ways to achieve the desired prescribed base 
motion. The dynamic force can be directly applied to the COSINE term of the 
ring No. 2 in Figure 6b. It can also be applied, using POINTAX cards, at 
selected points in the bottom layer. As observed in the test runs, both 
approaches gave the same results. The "POINTAX" approach was taken on the real 
model. 


Ike appropriate E/p ratio was determined to be of the order of 10 • This 
specific ratio gave an uniform base acceleration as desired* In addtition, 
both CONEAX and TRAPAX elements vere tested out zb the base layer. The TRAPAX 
element vas employed in the real model* 


Nodal transient analysis vas performed on this test model, t -ever, as 
noted from the eigen value analysis on the real model, there are a large number 
of natural modes in the low frequency range* Also, the impulse duration is so 
short that modes with extremely high frequencies are expected to participate in 
the dynamic response* Obviously, modal transient approach will not be accurate 
without using a large number of modes, which in turn will make it extremely 
time-consuming. Therefore, direct transient approach was taken for the 
analysis of real model* 


STATIC ANALYSIS 

Static analysis of both bonded and unbonded cases vere performed as checks 
on the fuze model* The load condition considered was that of 100 lbs. 
transversely applied at node 9. 

For this loading case* deflections at the load point and nodes 17, 21, are 
listed for different choices of harmonic terms* (See Table 2) Based on its 
convergence characteristics, AXIC = 3 was specified for the transient analysis* 
This means one constant and three harmonic functions were combined to describe 
the deformation field along the circumferential direction. Stress distribution 
in different parts vere examined, and the MPCAX cards were modifid accordingly 
to achieve appropriate constraints between nose, collar and plate. 

In order to calibrate the model, simple experiments were devised to 
perform static tests on the fuze nose cones. A force gage was used to apply a 
static load onto the tip of the fuze* The fuze was mounted in a rigid base and 
a dial indicator was used tc record the fuze deflection at the point of the 
load* The test setup is shown in Figure 7* 
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Loads were applied near the tip of the fuze, incrementing the load from 20 
lbs. to 140 lbs. and recording the deflection of each increment. This test vas 
performed for bonded sample of the fuze. Comparison of NASTRAN and experiment 
results is in Table 3. Note that the noticeable differences may have been 
caused by the fact that the tip of the force gage made a slight indentation in 
the surface of the nose cone. 


EIGEN VALUE ANALYSIS 

Eigen value analysis were performed to further check the appropriateness 
of the NASTRAN model. MPCAX condition; between nose-collar-plate were 
repeatedly examined and modified. The final arrangement of these conditions is 
listed in Table 4. This specific arrangement induces only relatively low 
stress in the plate while maintaining stability of the model. It is worth 
mentioning that in the seventh field of the MPCAX Card, the HID number should 
take he form of (Integer 20) instead of (Integer >0). Otherwise, undesirable 
rigid body modes may be present in the model. 

Ten fundamental frequencies and modes were extracted for both bonded and 
unbonded cases. As listed in Table 3, the first 9 eigen values of both cases 
are exactly the same. Also from Table 5, the lowest four frequencies are very 
close together at 30 cps. It was identifed, very interestingly, that each 
harmonic contributed to one of them. 

DIRECT TRANSIENT ANALYSIS 

Transient responses from t = 0 sec to t = 270 x 10 * sec, which was three 
times the l^nulse duration, were obtained for both bonded and unbonded cases* 
Solution convergence with respect to different time steps is listed in Table 6. 
at * 2 x 10‘* sec was chosen for the time increment in the analysis. The base 
displacement time history compared very well with the prescribed one. The 
dynamic stresses in the plane of impact in elements 104 and 103 were plotted 
and compared between both cases (see Figures 8-11). 

Figures 8 and 10 show the hoop stress time histories of both cases. The 
peak hoop stress in unbonded collar is up to about three times higher than that 
in the bonded collar. This explains the hoop failure in the unbonded collar as 
was observed from experiments. 

Figures 9 and 11 show the axial stress time histories of both cases. The 
peak axial stress in bonded collar is significantly higher than that in the 
unbonded collar. 

CONCLUSION 

The NASTRAN results presented in Figures 8 through 11 lead to the 
following conclusion: 

By bonding the nose and the collar, the nose joint will take dynamic load 
mainly in the form of axial stress. The peak hoop stress experienced in the 
bonded collar will thusly be reduced by a factor of 4. Hence, the hoop failure 
was effectively avoided in the bonded collar. 
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TABLE 1. SOLUTION CONVERGENCE WITH RESPECT TO MESH REFINEMENT 


MESH 

ASPECT 

RATIO 

TIP DEFLECTION 

BENDING 
STRESS 
AT POINT A 

8x2 

wm 

0.001832 

— 

16 x 2 

3.75 

0.003029 

798 

32 x 1 

0.94 

0.003568 

2151 

' 

! 

32 x 2 

' 

1.875 

0.003631 

2085 

64 x 1 

0.47 

0.003775 

2314 

ANALYTICAL 

— 

0.003918 

2292 
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TABLE 2. CONVERGENCE OF LATERAL DEFLECTION 


WITH RESPECT TO THE NUMBER OF HARMONICS (BONDED COLLAR) 


NO. OF 
HARMONICS 

LATERAL DEFLECTION 

NODE 9 

NODE 17 

NODE 21 

3 

-0.00676 

-0.00467 

-0.00375 

4 

-0.00737 

-0.00467 

-0.00375 

7 

-0.00767 

-0.00467 

-0.00375 
1 1 


TABLE 3. COMPARISON OF STATIC RESULTS 
(TRANSVERSE LOAD AT 0.27" FROM NOSE TIP) 


CASE 

EXPERIMENT 

NASTRAN 

BONDED 

-0.014 

-0.00676 

UNBONDED 

-0.020 

-0.00821 
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TABLE 4. CONSTRAINT CONDITIONS ON HARMONICS 
BETWEEN NOSE, COLLAR AND PLATE 


CONSTRAINED 

RING 

PAIRS 

BONDED 


UNBONDED 

r 

e 

z 

r 

9 

z 

35 

235 

0-3 

0-3 

0-3 

0-3 

0-3 

0-3 

37 

237 

0-3 

0-3 

0-3 

0-3 

0-3 

0-3 

1 

40 

240 

0-3 

0-3 

0-3 

0-3 

0-3 

0-3 

43 

243 

0-3 

0-3 

0-3 

0-3 

0-3 

0-3 

46 

246 

0-3 

0-3 

0-3 

0-3 

0-3 

0-3 

i-9 

249 

0 

1 

OJ 

0-3 

0-3 

0-3 

0-3 

0-3 

52 

252 

0-3 

0-3 

0-3 

0-3 

0-3 

0-3 

55 

255 

0-3 

C-3 

0-3 

0-3 

0-3 

0-3 

58 

258 

0-3 

0-3 

0-3 

0-3 


0-3 

60 

260 

0-3 

0-3 

0-3 

0-3 

i 

0-3 1 

63 

263 

0-3 

0-3 

' 

0-3 

0-3 

1 

I 0-3 : 

64 

364 

— 


0 

— 


i o 

65 

365 

0,1 

— 

0.1 

0,1 

1 

0,1 ; 

70 

370 

1 

— 

0,1 

— 

— 

0.1 
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EIGEN VALUE ANALYSIS BONDED NOSECONE 


REA'. EIGENVALUES 


NODE 

EXTRACTION 

EIGENVALUE 

IAN 

CYCLIC 

NO. 

CRDER 


FREQUENCY 

FREQUENCY 

1 

6 

3.6116205+04 

1.900A79E+02 

3.C247Q6 f+01 

“> 

5 

3.626*25E+0* 

1.904317E+02 

3.030815r+01 


7 

3.642207E+04 

1. 70S457E+0? 

3.037*026+01 

4 

* 

3.64B653E+04 

1.9101535+02 

3.040102E+01 

r 

2 

9. 19E3565+04 

3.03287EE+02 

^.3269765+01 

6 

1 

1.029579E+05 

3. 208706E+02 

5.106813E+01 

*7 

3 

1.2119225+05 

3.4312675+0? 

5.5406035+01 

* 

9 

2.20552°5+05 

4.6^63065+02 

7.474403E- 01 

Q 

10 

3. 234188^+05 

5 . 626992c +02 

9.05 1129E +01 

10 

9 

3.374114F+05 

5 . 808713E +02 

9.2*48555+01 


EIGEN 

VALUE ANALYSIS 

UN30NDED NOSEC ONE 






REAL E I G 

ENVALUES 

NODE 

EXTRACTION 

EIGENVALUE 

RADIAN 

CYCLIC 

NO. 

ORDER 


FREQUENCY 

FREQUENCY 

1 

6 

3.6118205+04 

1 .900479E +01 

3.0247065+01 

2 

5 

3.626^255+04 

1.9043176+02 

3.03081? :-+0 

3 

7 

3.6422076+04 

1.708457E+02 

3. 0374026+01 

4 

4 

3. 6486836+04 

1.910153E+02 

3.0401C25+01 

'j 

2 

9. 19B3 565+04 

3.0328785+02 

^.3269765+01 

6 

1 

1.029579E+05 

3.208706E+02 

5.1068135+01 

7 

3 

1.2119226+05 

3, ,512675+0? 

5. 5406035+01 

8 

9 

2.205529E+05 

* 6863065+02 

7 .474*035+0 1 

? 

10 

3.2341885+05 

5.6B6992E +02 

9.0511295+01 

10 

9 

3.246078E+O5 

5.6974375+02 

9 .0677525 + 01 


TABLE 5. FIRST TEN NATURAL FREQUENCIES OF 
FUZE EXTRACTED FROM NASTRAN 
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TABLE 6. CONVERGENCE OF TRANSTEN SOLUTIONS WITH 
RESPECT TO TIME STEP SIZE 


BASE MOTION 


At = 2 m sec 

at = 1 M sec 

At = 0.5M sec 

Acceleration At: 





1 t * 8 m sec 

6.7 x 10 5 

6.08 x 10 5 

5.98 x ID 5 

5.92 x 10 5 

t = IO u sec 


9.21 x 10 5 

9.13 x 10 5 

9.09 x 

Displacement At: 





t = 8 m sec 

1 

5.53 x 10 - 6 

3.77 x 10 ~ 6 

| 

3.47 x 1 C ~6 

3.4 x 10 ~ 6 

t = 10 u sec 

i 

— 

8.63 x 10 ~ 6 

i 

i 

8.17 x 10 " 6 

8.06 x 10 -6 
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Lgure4a - Finite Element Mesh 

Of The Fuze Assembly 
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Figure 5b- Axial Pressure Loading 
Condition 
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Figure 5c - Transverse Tip Loading 
Condition (8x2 Mesh) 
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Figure 6 - Test Model II 
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Figure 7. Test equipment setup 
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